Large- N c QCD and Pade Approximant Theory 
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^ . Abstract 

! 111 the large- N c limit of QCD, the expansion of the vacuum polarization at low 

energies determines the whole function at any arbitrarily large (but finite) energy. 
This result is an immediate consequence of the Theory of Pade Approximants to 
Qh! Stieltjes functions. 
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In the limit of a infinite number of colors, QCD simplifies and becomes a theory 
of stable noninteracting mesons 0. The number of these mesons is infinite because 
Green's functions have to reproduce the logarithmic dependence on energy present in 
the parton model. In fact, in spite of all the attempts made, the solution of QCD in the 
large N c limit has not yet been found. The large-iV c limit of QCD is still too difficult 
to be solved exactly. It seems to be interesting, therefore, to develop a systematic 
approximation to QCD in this limit. In this paper, we would like to point out that 
Pade Approximants are such an approximation, at least in certain cases. 

Given a function f(z), of complex variable z, with power series expansion about 
the origin given by 

N+M 

f(z) » > (!) 

n=0 

we shall define the Pade Approximant, R^(z), as the rational function 

R ^ Z ) = ( 2 ) 

satisfying the condition that its expansion about z = matches N + M + 1 terms in 
the expansion of the function f(z) in Eq. (JTJ). Notice that, without loss of generality, 
one can always choose the first coefficient bo = 1. 

In general, the convergence properties of Pade Approximants are much more difficult 
than those of, e.g., power series and it is not guaranteed that they will converge to the 
function f(z) in the combined limit N, M — > oo. This of course calls into question the 
reliability of the approximation and probably explains why they are not such a popular 
tool in Quantum Field Theory as compared to power series. 

However, there is a very important exception to this general situation. If the 
function f(z) is a Stieltjes function whose expansion about the origin has a finite 
radius of convergence R > , i.e. 

with 2/(t) a bounded, nondecreasing function 1 with finite real- valued moments given 
by 

fn = (-!)"/ r n du(r), (4) 
Jo 

then there is a theorem in the Theory of Pade Approximants which assures convergence 
for N = M + J, J > -1, as M -> oo in Eqs. (JT]0)j21. 

As we shall now see, the vacuum polarization function is an example of a Stieltjes 
function of the momentum variable. To simplify matters, let us take the isospin limit 
(m u = rrid) and define Ily(Q 2 ), with q 2 = —Q 2 , as 

- gV) n v(Q 2 ) = / d 4 x e l(}X (T{u llM d(x)d lu u(0)}) ■ (5) 



1 Notice that v{t) is not even required to be continuous. 
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As it is very well known, the function U V (Q 2 ) satisfies a once-subtracted dispersion 
relation given by 



O 2 r°° dt 

U V (Q 2 ) = U V (Q) - / — — ^ lmU v (t) , (6) 

where M 2 is the mass of the lowest vector state, the rho meson. In the large- iV c limit, 
which we have already taken in Eq. (jfi|L the resonance spectrum is made of infinitely 
narrow one-particle states. This is the reason why the lower limit is cutoff at M 2 in 
Eq. (JHJ). Of course, we are assuming here that there is a mass gap in large- N c QCD, 
i.e. M 2 p ± 0. 

For any finite N c , two-particle states also contribute, which gives the spectral func- 
tion a support all the way down to t = 4m 2 . However, the physics between the two 
scales m n and M p is already well understood in terms of Chiral Perturbation Theory 
It is the physics in the intermediate scales between M p and the onset of the QCD per- 
turbative continuum at a few GeV which remains difficult. In Chiral Perturbation 
Theory, all the physics at and above the first resonance mass is encoded in a set of 
low-energy constants which, in principle, depend on the dynamics of QCD, although 
the precise connection is unknown. As the papers in Ref. [HJ E] show, it is for this 
connection that the large-iV c expansion may become most useful. Therefore, in the 
following we shall only consider the N c — > oo limit. 

The basic observation is actually extremely simple. Making the change of variable 
t — 1/r in Eq. (JHJ), one can rewrite this combination as 



^2._ ny(0)-IV(Q 2 ) fH* dv{r) 
Wj " Q 2 ~Jo 1 + rQ 2 ' 



with du(r) = -Imny(r _1 )(ir with Imny(r _1 ) > 0, i.e. one can express the function 
$(Q 2 ) precisely in the form (jSJ), with R = M 2 p . The integral representation (J7|) assures 
that the combination $(Q 2 ) is a Stieltjes function of the complex variable Q 2 and 
allows one to make connection with Pade Approximant Theory. 

The expansion of $(Q 2 ) about the origin is related to the corresponding Chiral 
Perturbation Theory expansion 3\ and can be written like 

oo 

$(Q 2 )=]TL n (Q 2 r • (8) 

n=0 

In the chiral limit (m u ,mrf — > 0) the coefficients L n are given by the appropriate low 
energy constants which accompany the operators with chiral dimension 2n + 6 in the 
chiral Lagrangian. These coefficients are given in terms of the spectral function by the 
following integral 

Cm? 1 

L n = (-1)" / P drr n -ImlMr- 1 ) . (9) 
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Due to the upper limit in the integral representation (|7jL clearly the power series 
expansion (JBJ is convergent only for \Q 2 \ < M 2 We would like to emphasize that the 
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underlying reason for the existence of the simple power series in Q 2 given by (JBJ) below 
Mp, i.e. with no chiral logarithms, is the cutoff at M p at the lower end of the spectrum 
in (JHJ) which is in turn a consequence of the N c — > oo limit. As is well known, chiral 
logarithms are subleading effects at large N c [3J|nilH|. 

The N c —>■ oo limit of the vector spectral function becomes the highly singular 
distribution, i.e. 

oo 

-Imllv (t) = F nM 2 J(t - Ml) , (10) 

n=l 

with the index n labeling the increasingly heavier sequence of vector resonances, i.e. 
M n = M p , M p >, Mp", ... as n = 1,2,3,... Consequently, the combination (JJJ) is a mero- 
morphic function and can be written as the infinite sum 

*« 2 >-E^. ("J 

n=l ^ " 

in terms of the resonance masses, M n , and decay constants, F n . 

On the other hand, the Pade R^ +J (Q 2 ) can always be decomposed using partial 
fractions as 

M 

RT\Q 2 ) = P [J/M] (Q 2 ) + £ > (12) 

r=l ^ "I" 

where p[ J / A/ l(Q 2 ) is a polynomial of degree J when J > 0, and it vanishes if J < 0. The 
above mentioned theorem [2] states that the Pade R^ +J (Q 2 ) in Eq. (J2J) will converge 
to $(Q 2 ) as M — >• 00, for any ,7 > — 1. In fact, the convergence takes place for any 
finite value of Q 2 in the whole Q 2 complex plane except on the negative real axis 
—00 < Q 2 < —M 2 , where the poles of $(Q 2 ), i.e. the resonance masses, are located. If 
the Stieltjes function under consideration is multivalued (as is the case of the vacuum 
polarization function due to the presence of logarithms), the limit M — > 00 of the Pade 
sequence selects the branch which is real for Q 2 real and positiveflj. This is precisely 
the physical definition of $(Q 2 ). All the poles of R^ +J (Q 2 ) lie on the negative real 
axis with —00 < Q 2 < —M 2 p . Furthermore, they are simple poles and with positive 
residues jS]. 

It is quite remarkable that out of the low-energy chiral expansion given in ((SJ) 
one may reconstruct the full function $(Q 2 ). Notice, in particular, that this function 
contains logarithms of the variable Q 2 at euclidean high energies because the anomalous 
dimensions of the relevant operators in the operator product expansion do not vanish, 
in general, in the large- N c limit. The above convergence theorem means that all this 
complexity is reproduced by the rational approximant (|T2*|) as M ^ 00. Pades offer 
a proper mathematical way to approximate certain large-iV c QCD functions, such as 
the vacuum polarization $(Q 2 ) in ([lip, by means of a rational function, i.e. with a 
finite number of poles. We think that these results may give a new perspective to some 
phenomenological approximations used in the past, starting with the celebrated vector 
meson dominance |6j. 
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Some results about the rate of convergence are known. For instance, in the case of 
the R^~ l Pade one has that [I] 



S{Q*) - R»-\Q>) \<K 



, n rr , rr 2M 



(13) 



where K is an unspecified constant. The bound is valid for any complex Q 2 , provided 
it is not on the cut — oo < Q 2 < —M 2 . Notice that the numerator in the bound (fT3|) 
is smaller than the denominator and, therefore, it goes to zero as M — ► oo. 

A word of caution is required, however. The fact that the approximant ()12j) con- 
verges to the full function $(Q 2 ) for any finite value of Q 2 does not mean that one 
can expand this approximant in powers of 1/Q 2 at high energies. On the contrary, any 
such expansion cannot be exact as it is unable to reproduce the logQ 2 behavior due 
to the nonvanishing anomalous dimensions. In other words, the limits Q 2 — > oo and 
M — » oo do not commute, as it is clear from Eq. (|13j) . 

In this respect let us mention that a completely different subject is the use of the, 
so-called, two-point Pade Approximants 8J. These are also rational functions similar 
to (J2J), but for which the construction of the two polynomials is done by matching 
simultaneously both expansions of the function f(z) around z = and z = oo. In this 
case, if the function f(z) has a power series expansion in 1/z at infinity, these powers 
will of course be trivially reproduced by construction. However, as we have already 
emphasized, Green's functions in QCD have logarithms in the large Q 2 expansion. 
Strictly speaking, this presents an obstruction to the construction of these two-point 
Pades. In certain cases where the Green's function considered vanishes to all orders of 
perturbation theory 2 , these logarithms enter only as a correction to the leading 1/Q 2 
power and they are screened by at least a power of the coupling constant a s . It becomes 
then justifiable to neglect these logs in a first approximation jHj. In this case, however, 
the Pades no longer approximate the full function and their convergence properties are 
not known. 

In the case of the vacuum polarization function (|11|) . the fact that R^ +J (Q 2 ) con- 
verges when M — > oo means that the polynomial P^ J / M \Q 2 ) in Eq. (|12j) tends to zero 
in this limit. Consequently, out of all the possible Pade approximants, the R^~ 1 (Q 2 ) 
Pade is particularly interesting because it allows a direct comparison of its poles and 
residues with the resonance masses and decay constants appearing in the function 
$(Q 2 ) in Eq. ((TTJ) . 

Restricting ourselves to this Pade R^f 1 in Eq. (JT2J) (i.e. with pV/ M \Q 2 ) = 0), 
one can easily find relations among resonance masses and low-energy constants. For 
instance, one can obtain a first-order approximation to the p mass as the position of 
the only pole of the lowest-order Pade R!{{Q 2 ). This turns out to be located at 

»l = -Y- ( 14 ) 

This prediction can be refined by looking at the he next-to-leading Pade, R\(Q 2 ). Its 
two poles at /if and /x| yield approximate values for the p and p' masses which verify 

2 Such as, e.g., the left-right vacuum polarization function Ulr in the chiral limit. 
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Figure 1: (Left) Plot of the ratio 
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as a function of Q 2 , for M = 1,2,3 and M = 11. 



Notice the large values of Q 2 considered. (Right) Close-up of the same curves with M = 1,2,3 
in the low-Q 2 region. 



the relations 



l4 + f4 



L1L2 — LqLj 



1 2 



J_ _^ uiu 2 - ^0^3 r\ 

T T i|2 ,,2 T 2 T T • V i0 / 

The above prediction for the first two resonances is not equally good for both of them, 
however. Because the Pade approximation worsens as one moves away from the origin, 
one should expect that the prediction for a resonance mass is better the lower the 
resonance. This means better predictions for the p than for the p', and the latter 
better than for the p", etc... for a given Pade. 

Should more terms from the chiral expansion be known, we could in principle go 
on refining the above predictions in a systematic way. Regretfully, the L n are presently 
not known. Nevertheless, in principle these terms can be determined in the future by 
analyzing scattering processes of pions and kaons in Chiral Perturbation Theory or 
by other methods such as, e.g., the lattice. When this will be done, it will be interesting 
to check the above relations. 

In the mean time we shall content ourselves with a simple illustration of some of 
these results by means of a toy model. To this end, let us consider the following 
function defined by the sum[lOJ [TT] 



,2 ~"~ ,,2 



1 

t4 



LiL 2 — L L 3 



*To„(Q 2 ) 



E 



1 



•f n{Q 2 + n) ' ^ ^ 

as a very simple toy model for the once-subtracted vacuum polarization (JJJ) in the 
large-Nc limit. Obviously its spectrum consists of masses located at M 2 = 1,2,3,..., 
and one may use the mass of the lowest resonance as the unit of energy. An exact 
result for <&T oy {Q 2 ) can be readily obtained as 

l (Q 2 ) = ^\m 2 + v+i] (17) 



Toy{ 



Q 2 



where the function ip(Q 2 + 1) is the Digamma function, defined as 



d 

dz 



18) 
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4.3711 


2.2967 


2.0409 


2.0038 


2.0002 


M 2 


= 3 






10.1201 


4.1755 


3.2698 


3.0526 


Ml 


= 4 








19.5578 


6.9160 


4.8688 



Table 1: Squared masses of the resonances obtained as poles of the corresponding Pade. 
For comparison, the true squared masses are also shown on the leftmost column. 

with r(z) the Euler Gamma function and 7 ~ 0.577216 being the Euler-Mascheroni 
constant. The function &Toy(Q 2 ) has the "chiral" expansion 

00 

$to,(Q 2 )^^(-1) p C(p + 2)Q 2p , (19) 

p=0 

where ((p + 2) = Y^Li n ~ p ~ 2 is the Riemann's ( function. On the other hand, the 
large Q 2 expansion is given by (| argQ 2 | < 7r) 

<fTo,(g) ^ +^ 2^ 2p ' (20) 

where 5 2p = (— l) p+1 C(2p) are the Bernoulli numbers for p = 1, 2, 3, .... 

While the chiral expansion has a finite radius of convergence \Q 2 \ < 1, the large 
Q 2 expansion has zero radius of convergence (i.e. it is asymptotic) due to the factorial 
growth of the Bernoulli numbers. Notice also the presence of the logQ 2 behavior in 
(J2UJ). Even so, Fig. 1 shows how the Pade Approximants -R*/" 1 are able to reproduce 
the full function at any finite value of Q 2 . Notice that no step function simulating a 
perturbative "continuum" is needed jllj . The fact that they approach the true function 
from below is also well-known j^j. 

As to the prediction for the resonance masses, Table 1 shows the position of the 
poles of the given Pades and compares them to the true values for the masses. As 
anticipated before, the prediction for the masses gets worse the heavier the resonance, 
for a given Pade. The prediction for a low resonance, however, is certainly reasonable 
even for a low-order Pade and becomes much better as the order of the Pade increases. 

In conclusion, Pade Approximants, when properly used, may turn out to be very 
useful tools for studying large-iV c QCD. We have exemplified this in the case of the 
vacuum polarization because, having a positive spectral function, it happens to be a 
Stieltjes function, to which a known convergence theorem applies. But this case is 
not unique. Similar results will also apply to other correlators with positive spectral 
functions. Apart from this, it would certainly be most interesting to find out under 
which conditions Pade Approximants converge in more general situations when the 
spectral function is not positive definite. 

The so-called two-point Pades also seem to be very promising, as they are able to 
connect the Operator Product and Chiral expansions j^]. However, to our knowledge, 
the consequences of their convergence properties for a QCD Green's function have not 
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yet been analyzed. It is clearly tempting to speculate about the possibility that they 
may be linked to theories with extra dimensions |12|. 
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